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Abstract 

We describe in detail the constrained procedure of differential renormalization 
and develop the techniques required for one-loop calculations. As an illustration we 
renormalize Scalar QED and show that the two-, three- and four-point Ward identities 
are automatically satisfied. 



1 Introduction 

Differential renormalization (DR) [Q] is a method of renormalization in coordinate space 
that yields directly finite Green functions without intermediate regularization or explicit 
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counterterms. It has proved to be quite simple and powerful in a number of applica- 
tions Q (see also Refs. |5| for formal developments). Standard DR manipulates 
singular objects as if they were well-defined, expresses them in terms of simple singular 
functions, and substitutes these by their renormalized value. The renormalized func- 
tions contain arbitrary dimensionful constants which play the role of renormalization 
scales and carry all the ambiguity inherent to the formal manipulation of singular ex- 
pressions. In symmetric theories, these scales can be adjusted in such a way that the 
Ward identities for renormalized Green functions are satisfied. For practical purposes, 
however, one would prefer a procedure that directly rendered symmetric expressions, 
without having to impose the Ward identities after each calculation. 

Such a procedure, constrained differential renormalization (CDR), was introduced 
in Ref. Q at the one-loop level. The idea is to impose that the renormalized expres- 
sions be compatible with a minimal set of consistent formal manipulations (rules). As 
a result, the ambiguities and arbitrary renormalization scales of DR are fixed and the 
resulting renormalized Green functions automatically preserve Ward identities. This 
was explicitly shown in Ref. Q for a variety of cases in abelian gauge theories. Su- 
persymmetry was also preserved in the calculation of supergravity corrections to the 
anomalous magnetic moment of a charged lepton |0, |8|. The method also works for non- 
abelian gauge theories, as shown in Ref. |^]. The extension of the constrained method 
to higher orders is more involved and needs further study. At any rate, it requires a 
good understanding and systematization of the one-loop order. 

The purpose of this paper is to describe CDR in detail and introduce the techniques 
required for one-loop calculations in any renormalizable theory in four dimensions. As 
an example, we renormalize (to one loop) Scalar QED, which contains all possible kinds 
of diagrams and derivative structures that appear in renormalizable theories. The plan 
of the paper is the following. Section 2 contains the procedure and rules of CDR. We 
first motivate them taking as a guiding example the vertex Ward identity in QED. 
Then we discuss how to expand any one-loop Feynman graph in a complete set of 
basic functions and we give the renormalization rules. Finally, the rules are used to 
determine the renormalization of the singular basic functions. The reader interested in 
practical calculations only needs the tables in this section. The full derivation is given 
in Appendix A. In Section 3 we apply the method to Scalar QED. We renormalize the 
two-, three- and four-point IPI Green functions and verify that they satisfy the Ward 
identities. Section 4 is devoted to conclusions. In Appendix B we collect the Fourier 
transforms of CDR expressions. An efficient computer code performing all operations 
automatically (in momentum space) is available |1C]. 



2 Constrained differential renormalization 

In general, the bare expressions of Feynman diagrams in coordinate space are too 
singular at coincident points to behave as tempered distributions. The role of renor- 
malization is to consistently replace these "singular" expressions by "regular" ones. To 
carry out this program, we first reduce each Feynman graph to a sum of (singular) 
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basic funtions and their derivatives, and then define the renormaUzed diagram as the 
corresponding sum of the renormahzed basic functions. CDR is essentially a set of rules 
for manipulating the singular expressions and fixing the renormalization of a complete 
set of basic functions. These rules are chosen in such a way that Ward identities are 
respected. We specify the rules below, after motivating them with a simple example. 
Then, we define the (one-loop) basic functions and apply the rules to renormalize them. 



2.1 Motivation of the rules 

Ward identities among Green functions are derived from the symmetries of the action 
using general properties of Quantum Field Theory, as gathered in Schwinger's action 



principle [11|. These properties are only formal for bare Green functions because they 
involve ill-defined expressions. Regularization and renormalization can invalidate them, 
breaking the symmetries at the quantum level. In absence of anomalies, the symmetry 
can be restored adding adequate finite local counterterms, as a consequence of the 
quantum action principle |12, However, in real calculations it is preferable to have 



an invariant procedure rendering directly symmetric expressions. This is the case of 
dimensional regularization with minimal substraction in gauge theories 15 1. In 



perturbation theory, the fulfilment of the Ward identities (and of the action principle, 
in general) relies heavily on the following property |15, the application of the kinetic 



differential operator to the propagator corresponding to some line in a Feynman graph 
is equivalent to the contraction of the line to a point. In other words, the fact that the 
free propagators are Green functions of the corresponding quadratic lagrangian must 
not be spoiled by the renormalization procedure. Actually, the corresponding proof 



for gauge theories in Ref. |15] is based on the fact that this property holds for the 
dimensionally regulated diagrams. We shall see that the propagator equations also 
play a central role in CDR. Although no explicit regulator is available in DR, one can 
impose simple conditions on the renormahzed expressions that guarantee the validity of 
the propagator equations and of other formal identities. To see what kind of constraints 
one must impose, let us consider the case of the vertex Ward identity in massless QED, 

(5^ - dl)V^.{x, y) = ie[S(x)5(y) - S(y)5(x)] , (2.1) 

where is the IPI photon-electron-electron vertex, S is the electron self-energy and 
X = xi — X2, y = X2 — X3, with 2:1^2,3 the external points (x2 being the point attached 
to the photon)Q For bare functions this Ward identity is easy to derive formally at 
one loop. Using the Feynman rules in Ref. |l^ (we work in euclidean space and in the 



^For convenience when treating the general case, we use in this paper the variable y with a minus 
sign with respect to the one in Refs. ^ ^, p^ . 
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Feynman gauge), 

(9^ - dy)V^{x,y) = {dl - ap[(-ie)37,r A(x)7^^^A(y)7„A(x + y)] 
= ie^[7„ r A(x) ^s^A(y)7„ A(x + y) 
-7„r A(x) ^s^A(y)7„A(x + y) 

+7ar A(x)7^r A(y)7,(5^ - a|l)A(x + y)] (2.2) 
= ie3[-J(x)7„^2/A(y)7^A(x + y) 

+%)7«rA(x)7„A(x + y)] 
= ieM(x)E(y)+<5(y)S(x)], 

where A(x) = is the massless Feynman propagator, stands for and 

□ = d^d^. The manipulations used in Eq. ( |2.2| ) are only formal because we are deal- 
ing with singular expressions (at x = y). We have first performed the linear change 
of variables described above, and used the Leibnitz rule to make the external deriva- 
tives act on individual propagators. Then we have used the propagator equation for 
a massless fermion (which is equivalent to replacing by □ and using the scalar 
propagator equation, DA = —6). The symmetry of A(x + y) under the interchange 
of X and y makes (5^ — 5pA(x + y) vanish; finally, the definition of the Dirac delta 
function justifies the last equality. When these ill-defined expressions are regularized 
or replaced by finite ones (renormalized) , some of these manipulations may not be jus- 
tified. For example, if the fermion propagator were naively regulated a la Pauli-Villars, 
5^'^(x) = ^^(A(x) - Aa(x)), with Aa the Feynman propagator of mass A, the prop- 
agator equation would be modified: ^S^^ = — A^Aa(x), and the third equality in 
Eq. ( |2.2| ) would no longer hold. Other two examples studied in the DR literature are 
the regulators introduced in Refs. @ and respectively. In the former, the modi- 
fication of the propagator invalidates the propagator equation; in the latter the short 
distance expansion invalidates the change of variables and the direct use of the delta 
function^. Our aim is to find a DR scheme automatically preserving the Ward identities 
for renormalized amplitudes. One may require the whole Eq. (|2.l| ) to be satisfied by the 
corresponding renormalized expressions, as we essentially did in Ref. Q to relate dia- 
grams with different topology. However this has to be improved because diagrams are 
expected to be renormalized independently of the particular combination in which they 
appear in the Ward identities. Instead, we demand that the manipulations performed 
to derive Eq. ( |2.2D be always valid for renormalized expressions. In other words, we 
require that DR commutes with this kind of manipulations. In what follows we show 
how this can be done consistently at one loop. Indeed, a few rules will be sufficient to 
fix the ambiguities of standard DR and to ensure the validity of the Ward identities at 
the same time. 

2.2 Rules for constrained differential renormalization 

We distinguish two kinds of rules. The first ones simply state that renormalization 
commutes with algebraic identities that allow to perform "straightforward" manipu- 

^ Both methods were adequate to study the relation between DR and the usual counterterm ap- 
proach, although the spirit of DR is precisely to avoid the use of such regulators. 
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lations, like sums of terms, Dirac algebra (in fom' dimensions), or application of the 
Leibnitz rule. With these manipulations one can express any one-loop Feynman dia- 
gram in terms of a set of basic functions that will be defined in the next subsection. The 
rules of the second kind extend well-defined identities of distribution theory to more 
singular expressions. They will be used to determine the renormalization of the basic 
functions. Let us enumerate them (the first two rules are essentially the prescriptions 
of the method of DR) [|: 

1. Differential reduction: singular expressions are substituted by derivatives of reg- 
ular ones. This can be done in two steps: 

(a) Functions with singular behaviour worse than logarithmic ('^ x~'^) are re- 
duced to derivatives of logarithmically singular functions without introduc- 
ing extra dimensionful constants. 

(b) Logarithmically singular functions are written as derivatives of regular func- 
tions. The usual DR identity p| 



1 logx^M 
4 X 



□ , (2.3) 



is sufficient to one loop. This identity also applies to massive expressions 
when they are expanded in the mass. Eq. (2.3) introduces the unique di- 



mensionful constant of the whole process, M, which has dimensions of mass 
and plays the role of the renormalization group scale. 

2. Formal integration by parts: derivatives act formally by parts on test functions. 
In particular, 

[dFf = dF^ , (2.4) 
where F is an arbitrary function and R stands for renormalized. 

3. Delta function renormalization rule: 
[F{x, XI, ...,xn)d{x - y)f = [F{x, xi, ...,Xn)]^S{x - y) . (2.5) 

4. The general validity of the propagator equation: 
F(x,a;i,...,x„)(n^ -m2)A„(a;)] = [F(x, xi, Xn)(-5(a;))]^ , (2.6) 

where Am{x) = mKi{mx) ^ modified Bessel function 

Rule 0, combined with rule reduces the "degree of singularity" , connecting singular 
and regular expressions. Forbidding the introduction of dimensionful scales outside 
logarithms, we completely fix the scheme. The last three rules are valid mathemat- 
ical identities among tempered distributions when applied to a well-behaved enough 
function F. The rules formally extend their range of applicability to arbitrary func- 
tions. They require that these identities commute with the process of renormalization. 
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Rule |2| is essential to make sense of rule for otherwise the right-hand-side of it would 
not be a well-defined distribution. Rule ^ connects functions with different number 
of propagators and has actually been applied in the literature to reduce three-point 
functions to two-point functions, thus allowing the use of the DR identity (^.j). As we 
shall see in the next subsection, it has further important implications. Note that the 
massive propagator equation ( |2.6| ) carries the necessary information for the propagators 
appearing in usual theories. Indeed, the corresponding equation for a massless scalar 
propagator is just the limit m — > of Eq. (^), the equation for a fermionic propagator 
follows from it, after the use of = □, and the equation for the propagator of a gauge 
boson in a general covariant gauge can be derived from the first two terms of the mass 
expansion of Eq. (|2.6|) (see Ref . ) . 

All these rules allow to renormalize any one-loop Feynman graph. Other possible 
manipulations can be incompatible with them and introduce ambiguities. One can still 
perform them as long as one keeps track of the arbitrary local terms they introduce. 
These local terms can then be fixed imposing consistency with the renormalization 
rules. In particular, two dangerous operations are explicit differentiation (which is 
not well-defined at coincident points)^ and index contraction ||6|, [l^, as we shall see 
explicitly later on. These two manipulations do not commute in general with CDR. 

The actual procedure of renormalization involves two steps: (1) express a Feyn- 
man diagram in terms of basic functions and (2) replace the basic functions by their 
renormalized value. The renormalization of the basic functions is done once and for 
all and does not depend on the particular calculation. This is the subject of the next 
subsection. Let us now comment on the first step. At one loop, a Feynman graph in 
coordinate space contains products of propagators with differential operators acting on 
some of them and, possibly, delta functions. It can also include constant objects like 
group factors, gamma matrices, the metric and the Levi-Civita tensor. To evaluate 
the diagram, we first make a convenient change of variables (see below). Then, we 
perform all the (Dirac) algebra. In particular, we do all possible contractions. This is 
a necessary prescription because index contraction does not commute with CDR. This 
prescription leads to universal renormalized basic functions. Finally, the Leibnitz rule 
is used to reorder the derivatives in each term, so that they act either on the whole 
product or on the last propagator. Symmetries among the space-time points can also 
be safely used as, for example, in F{x, y)dj{A{x + y) = F{x, y)d^A{x + y). 

2.3 Renormalization of basic functions 

We shall assume that a Feynman-like gauge is chosen for the gauge fields, so that their 
propagators are proportional to the scalar Feynman propagator. The formalism can 
be directly extended to general covariant gauges, as was indicated in Ref. Up to 
delta functions, any one-loop IPI graph is a linear combination of products of scalar 
propagators Ami{xi — Xj+i), with differential operators O^' acting on them. All points 
Xi are external and appear in a cyclic way. Therefore, using the Leibnitz rule, each 

^ Taking derivatives is essential to obtain the DR identities, but the corresponding ambiguities are 
eventually fixed by the CDR rules. 
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term of the graph can be written (ignoring constants and delta functions) as a sum of 
total derivatives of the whole product of propagators, with all the internal derivatives 
acting on only one adequately chosen propagator: 

O^^Am.ixi - X2)0^'A„2(X2 - Xg) • • • 0^"A^Jx„ - Xi) 

:axt{i)[Ami(2i)A„,(z2) • • • Of^,(.)A„„(zi + Z2 + . . . + . (2.7) 



In the first line we have performed a convenient change of variables: zi = xi — X2, 
Z2 = X2 — X3, . . . , Zn-i = Xn-i — Xn, Zn = Xn- This eliminates one variable (zn) due to 
translational invariance. Besides, since every Zi appears only in one of the first n — 1 
propagators and in the last one, it is straightforward to use the Leibnitz rule to make 
all derivatives act on the last propagator and obtain the second equality. Because of 
rule ^, the renormalization of the graph reduces to renormalizing expressions of the 
form 

A^,(zi)A„2(22) • • • A„„_,(z„_i)0^i A„„(zi + Z2 + --- + zn-i) , (2.8) 

which we call basic functions. For massless basic functions we shall suppress the mass 
subindices. In renormalizable theories, singular one- loop diagrams involve basic func- 
tions with at most four propagators and no more derivatives than propagators. All the 
singular basic functions for this class of theories (in the Feynman gauge) are displayed 
in Table 1. We denote A, B, T and Q the basic functions with one, two, three and four 
propagators, respectively. In the following we shall use x, y, z to denote 2:1,2,3, respec- 
tively, and assume that they are the (ordered) arguments of the basic functions unless 
otherwise specified (e.g., T[n] = T[n](x,y)). The renormalization of these functions 
is carried out using the rules of CDR Renormalized one- and two-point functions 
are gathered in Tables 2 and 3, and three- and four-point functions, in Table 4. For 
the reader's convenience the massless A and B functions are given in Table 2, although 
they can be obtained from the massive ones in Table 3, taking the appropriate limit. 
The expressions for T and Q functions in Table 4 apply directly for any value of the 
masses. The Fourier transforms are collected in Appendix B. The reader interested in 
applications only needs the expressions in Tables 2, 3 and 4 for coordinate space or 
Tables 5, 6 and 7 in Appendix B for momentum space calculations. 

Note that it is important in our procedure to distinguish basic functions with con- 
tracted and uncontracted differential operators, because contraction of Lorentz indices 
does not in general commute with CDR. For instance, from Table 4 |^, 17 1, 



T^[n] = [S^,T[d^d,]f + b^,T''\d^d,\ . (2.9) 

This justifies our prescription of contracting indices before identifying the basic func- 
tions. In the following we briefly describe how these tables have been obtained. The 
detailed derivation can be found in Appendix 
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logarithmic 


linear 


quadratic 


cubic 


1 prop. 






A^[l] 




2 props. 




Bmim2 [^m] 


Bmim2 [°] 
Bmim2 [^^^J^] 




3 props. 


Tmim2m3 ['— '] 
Tmim2m3 ["9/^5;/] 








4 props. 


Qm,im2m3rri4 ['—"—'] 
Qmim2m3m4 'f^/^f^;/] 
Qmim2m3rri4 [c^/i^^i/f^pf^cr] 









Table 1: Singular basic functions for renormalizable theories in four dimensions. Lines 
are ordered according to the number of propagators and columns according to the 
degree of singularity. The function A, that appears in tadpoles, is defined as Am[C] = 
OAm{x)S{x). 

In general the renormalization of basic functions proceeds in two steps. First, a 
differential equation is solved for non-singular points, in order to express the singular 
basic functions as derivatives of well-behaved functions. A useful trick in the case of 
complex tensor structures is to decompose them into trace and traceless parts: the 
former carries the leading singularity but is simpler, while the latter is less singular. 
Second, the arbitrary local terms are determined according to the CDR rules. Rule || 
is actually just an initial condition for the local terms, while the remaining rules relate 
the local terms of different basic functions. Basic functions with simple enough tensor 
structures can be directly expressed in terms of other basic functions that have been 
previously renormalized. 

Let us consider first the renormalization of Bm,im2[l] aiid Am[l]. We start with 
the massless case. Direct application of rule |l] gives 



B^[l] 



(47r 



2^2 



^ 1 logx^M^ 



which is the standard DR identity. The fact that A[l] = A{x)6{x) is local, together 
with power counting, implies that the most general renormalized value of this function 
is of the form 

A^[l] = (an + ^i^)^^^)^ (2.11) 



Rule la tells us not to introduce the dimensionful constant /i, so the second term in 



the equation above vanishes. Eq. ( 2.1C| ) and = in Eq. (|2.11| ), based on rule ||, are 



the initial conditions of the renormalization process. This rule is not needed anymore. 
As it is shown in Appendix rule ^ implies that the first term in Eq. ( |2.11| ) must also 
vanish, so we have 

A^[1]=0. (2.12) 
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A^[l] = 

A«[a^] = 

B^[l] = -6^°^ 
B^[n] = 



Table 2: Renormalized expressions of massless one- and two-point basic functions. 

For massive basic functions we use recurrence relations among modified Bessel functions 
to obtain the expressions for non-singular points (see Ref. ||5|, |l8| and Appendix C of 
Ref. §). The local terms of B^^^Jl] and A^[l] are fixed so that in the massless limit 
we recover B^[l] and A-^[l], respectively. In this way we find: 



R 



-Bjni mo i'^] 



1 [ miKi{mix)m2Ki{m2x) 



(47r2)2 

1 f 17111712 ^ .2T^oi'mix)Ki{m2x) + Ko{m2x)Ki{mix) 

'□ - (mi -I- 1712) \- 



32-71^ \ nil + "1-2 X 
t2 



+ 2vr^llog^ + !^log^U(.)^ (2.13) 
\ 77111712 mi + 1712 "mi J J 

where M = 2M/^e and 7^ = 1.781 ... is Euler's constant. The massive one-point 
function A,^[l] is determined (see Appendix ^) from A^[l] and B^^^^: 

1 /If 2 

A^[l] = -^rn\l - log -^)6{x) . (2.14) 

The renormalization of the remaining basic functions is obtained from B^^^^ll] and 
Am[l]^ (B^[l] and A^[l] in the massless case) by recurrence relations based on rules ^, 
^ and 0. These recurrence relations follow from the Leibnitz rule and two operations: 
point separation and point contraction. Point separation ||^] allows to relate a generic 
renormalized basic function with n propagators and r derivatives, Fmi^.m„ [OW], to 
renormalized basic functions with n + 1 propagators and r, r -f 1 and r + 2 derivatives. 
Using the rules, 

¥t\'^.mM{zi,...,Zn-l)5{Zn) 
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167r2 

= 



nR [11- if mima / , \2i Koimix)Ki(m2x)+Ko(m2x)Ki(mix) 



+ 27r2 f log + ^^^1^^ log ^) S(x) ] 



m2-R"i(m23:) 

X 



mimal J 



"^i^^mima ^mi [1] 



+ 3(9;.5,-iV°)B^H™2[l] 



'rniKi(mix) 



+ 



1 

T67F 



} 



+ mlKo{m2x){d^di, - jd^^O) 



Table 3: Renormalized expressions of massive one- and two-point basic functions. 



[Amjizi) ... A 

[Am^Zi) . . . Am„_j(z„„i)0'=iAm„(zi H h Zn-l + Zn)S{Zn)]^ 

R 

- A^^ (zi) . . . A^„_j (z„_i)0^i A^„(zi + • • • + z„_i + z„)(D^" - m^+i)A^„^, 



(2.15) 



On the other hand, point contraction relates a renormalized basic function with a 
d'alambcrtian and r derivatives to renormalized basic functions with one less propagator 
and 0, 1, . . . , r derivatives: 



R 



\Z\_ ) ■ ■ ■ ^run [Zfi 



+ Zr, 



n+l^nin+i {Zl + ■ ■ ■ + Z; 



■n)) 



R 
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Tmim,2m3 [^M^i^] ~ Tmim2m3 [9/j<9,y — ^(^^^^D] + J^^i/T^^^^ms ['-'] ~ 1287r^ ^t^u^{x)S{y) 

T^.m^msPdf^] = mlT^.m^mAd,] " K^m^i^M^ + v) ' df, (b«^^J1]<^(x + y)) 

'^mim2m3 [^M^i^^p] ~ T^mim2ms[dp,d^dp — ■^{5pi,dp + Sp^pd^ + 5;yp9^)n] 

Qmim2m3m4 ['-"-'] ~ fn\^mim2m.-j,m.i ~ (Tr)iim2m3 [Ijl^ls^ + y + ^)) 
mim2m3 

Qmim2m3m4 [^^M^*^] ~ ''^4Qmim2m3m4 [c^/^C^!/] ~ 9^9^ (Tmim2m3 [ll^^C^; + J/ + 2^)) 

[d,]6{x + y + z))-dt (T 

mim2m3 

[dp]5{x + y + z)) 

- '^m^m2mz[dpdu\5{x + y + z) 
Qmim2m3m4 [^M^J^^P^cr] — ^m\m2imzva4\dpdi,dpd(j 24 ('^A^^'^P"" ~^ ^pp^ua ~\~ ^ pa^vp)^^ 

+ M^pu^pa + SppS^^ + 5paS^p) (Q^imamgmjnn] + ^d{x)6{y)6{z)'^ 



Table 4: Renormalized expressions of three- and four-point basic functions. The basic 
functions with traceless differential operators are directly finite. 



= ^ Ot {Ftl''.m,Mi^i, zn^,)5{zr + • . . + z„)) 

i 

+ [0]{zi,. ..,zn). (2.16) 

Oi (OQ is a differential operator with (r — rj) derivatives, < < r. Observe that 
point contraction and point separation are not inverse operations of each other due to 
the different delta functions in the l.h.s. of Eq. ( p. 15 ) and in the r.h.s. of Eq. ( 2.1(i| ). 



This allows to obtain non-trivial information when Eq. ( |2.16| ) is inserted in the last 
term of Eq. ( |2.15| ). In Appendix ^ we show how the combined use of the operations 
just described fixes the renormalization of all the singular basic functions of Table 1. 

3 Renormalization of Scalar QED 

In this section we apply the procedure of CDR to Scalar QED and verify that the Ward 
identities among renormalized Green functions are satisfied. The presence of derivative 
couplings makes Scalar QED the simplest renormalizable theory that contains all the 
singular basic functions in Table 1. In this context, the techniques described above can 
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Figure 1: Feynman rules of Scalar QED. Wavy lines correspond to photons and solid 

lines to scalars. In the trilinear coupling, (?/i=9^ — 9/*, with dn (dfj.) acting on the 
incoming (outgoing) scalar. 

be fully illustrated. We calculate to one loop the IPI Green functions of two, three 
and four points, which arc the only ones that require renormalization. The one-point 
Green functions directly vanish. In the following, after giving the lagrangian and the 
Feynman rules, we briefly discuss the renormalization of the six IPI Green functions 
that contain singular graphs, namely the vacuum polarization, the scalar selfenergy, 
the photon-scalar-scalar interaction vertex, and the three kinds of four particle vertices: 
photon-photon, photon-scalar and scalar-scalar scattering, and the corresponding Ward 
identities. We also check that the one-loop /3-functions and anomalous dimensions of 
Scalar QED are recovered. 

The complete lagrangian for Scalar QED in euclidean space and in the Feynman 
gauge is 

Ce = \f^uF^, + ^{d^^A^f + {d^ - ieA^)(P\d^, + ieA^)(l) + m^V + ji^pUf ■ (3-1) 
The corresponding Feynman rules are depicted in Fig.l. 
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Figure 2: Fcynman diagrams contributing to the vacuum polarization (a), scalar self- 
energy (b) and photon-scalar-scalar vertex (c). 
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3.1 Vacuum polarization 

The two Feynman graphs contributing to the vacuum polarization, al and a2, are 
depicted in Fig. 2. Their contributions are 



^(b)M = -2eH^,AUx)5{x) , (3.2) 



where x = xi — X2 is the coordinate difference. In the following, unless otherwise stated, 
we shall use the shifted variables defined in the text. In terms of basic functions the 
two contributions read 

^{b)tiiy{x) = -2e^5^ykra ■ (3.3) 

Using the expressions of the renormalized functions in Table 3, the sum is ]|6| 

„ ,- _N [1/ , 2\ /^-^o("^2:)ifi(mx) 



(□ - W)(- 





1 4 

+ rn^iKlimx) - Kf{mx))) + -7r2(log - -)5ix) 

3 3 



(3.4) 



which is transverse, as required by the Ward identity 



dm^,{x) = 0. (3.5) 



The scale-dependent part is 

M^n^,(x) = {d^d, - 'J(^) • (3-6) 

3.2 Scalar selfenergy 

The one^oop scalar selfenergy is obtained from graphs bl, b2 and b3 in Fig. 2: 

S(,)(x) = -e^ [□ (A(x)A„(x)) + 2do, (A(x)a„A™(x)) 

+A(x)nA„(x)] , (3.7) 

S(fe)(x) = -4e2A(x)5(x), (3.8) 

S(e)(x) = -AA^(x)5(x), (3.9) 

where derivatives acting on the external scalars have been changed by minus deriva- 
tives acting on the whole amputated expression. In terms of basic functions, the total 
renormalized selfenergy is 

Sfi(x) = -e2 (QBlil] + 2dpBUdp] + B^a] + 4A^[1]) - AA^[1] . (3.10) 

The explicit expression can be obtained using Tables 2 and 3. The scale-dependent 
part reads 
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3.3 Photon-scalar-scalar vertex 

The contributing diagrams, cl-c4, are shown in Fig. 2. Their renormahzed expressions 
in terms of basic functions are 

+ A{(r^ + di)d^-dym} , (3.12) 

V(f)^ = -2ie^ [dlBUl]{y) + Bo«„[5^](y)} 5{x) , (3.13) 

^(5^ = -2ie3 {a^Bl[l](x) + Bl[9^](x)} 5{y) , (3.14) 

^(5/.=0- (3-15) 

Graph c4, which is proportional to A, vanishes directly because of charge conjugation. 
The renormahzed vertex results from the sum ofEqs. (|3T^|3T^ ) and the use of Tables 3 
and 4. The scale-dependent part is 

The renormahzed vertex function and selfenergy fulfil the Ward identity 

{dl - a|l)K«(x, y) = ie{T.\x)6(y) - S(y)«5(x)) , (3.17) 



analogous to Eq. (2.1) in QED. This can be checked integrating with any test function 
(p{x,y). In particular, (j){x,y) = e*^'^ leads to simple expressions for both sides of the 
identity |jl^. Observe that the part of the Ward identity proportional to A (coming 
from diagrams b3 and c4) is trivially satisfied, for both sides vanish. In general, the 
diagrams proportional to A form a gauge invariant subset. 

3.4 Four-point functions 

The diagrams contributing to the photon-scalar (V^i/), dl-d7, and photon-photon (V^^pa-), 
el-e3, and to the scalar-scalar (V), fl-f9, scattering are depicted in Figs. 3 and 4, respec- 
tively. As in previous cases, the renormahzed Green functions are obtained expanding 
in basic functions and substituting these by their renormahzed expressions in Tables 2- 
4. The final results are too lengthy to be reported here. We have calculated them 
with a symbolic program based on Mathematica. The scale-dependent parts are readily 
obtained. For the photon-photon-scalar-scalar vertex we have 

4 



M^<,(x, y, z) = 5^,-^5{x)5{y)5{z) . (3.18) 



Notice that there is no A-dependent part. Actually, Eqs. ( 3.16| ) and ( |3.18| ) follow from 



gauge invariance and Eq. ( |3.6| ). The scale-dependent part of the four-photon vertex 
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cancels when all diagams are summed, as required by renormalizability. Finally, the 
scale-dependent part of the four-scalar vertex is 

M^y(x, y, z) = Y^^e^ - ^e^X + 5X^)6{x)6{y)6{z). (3.19) 

Observe that there is a A-independent piece, reflecting the well-known fact that a quartic 
scalar self-coupling is needed in the lagrangian for multiplicative renormalizability. The 
Ward identities for the four-point functions are more transparent in terms of the original 
external points: 

d1i^V^ixi,X2,X3,X4,) = ieVj^{xi,X4,X2){5{xi - X3) - 6{x2 - X3)) , (3.20) 
d^^W^,pAxi,X2,xs,Xi) = 0, (3.21) 

and the identities related to Eqs. ( |3.20| ) and ( |3.2l| ) by Bose symmetry. There is no 
Ward identity relating V to other singular Green functions. In terms of the shifted 
variables, Eqs. j ^lUj ) and (|3]2|) read 

id; - d^W^^Ax, y, z) = ieY^ix + y + z-y- z)[5{x + y) - 5{y)] , (3.22) 
5^^,Sp<x(^,y,^) = 0. (3.23) 



The scale-dependent parts of VJ^, and V^p^ trivially satisfy Eqs. ( ^.221) and (|3.23| ). 



The merit of CDR is to give also the correct scale-independent parts of the Ward 
identities. We have checked numerically that this is indeed the case, supplementing 
the program used for the calculation of the renormalized Green functions with the FF 
package |^9| {FF works in dimensional regularization, but we have only used it for finite 
functions). 

3.5 Re normalization group 

To conclude this section we compute the first order anomalous dimensions and (3- 
functions. The renormalization group equation for any renormalized amplitude reads 

^7^7 + f^--^ + f^^-^ + Imm^ - - ua^a) r^("^'"^) = , (3.24) 
oM oe OA dm J 

where j^^a are the anomalous dimensions of the corresponding fields and F-^^""*'"^) is 
the IPX renormalized amplitude of scalar and ua photon fields^. The tree level part 
of the F functions can be read off from the Feynman rules in Fig. 1 and the scale- 
dependence of the one-loop amplitudes is given in E qs. (PI) , ( pT| ), ( |3l^ ), ( |3l8|) and 



(|3.19| ). Using the renormalization group equations (|3.24| ) for the relevant amplitudes 
one obtains 

= (3.25) 



* The only effect of the gauge parameter dependence is to consistently cancel the tree-level longitu- 
dinal term in the equation for the vacuum polarization. 
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e2 



7</. = -^, (3.26) 



7r2 

2 



= (3.27) 
A = (3.28) 



24e^ - 12e^A + SA^ 
167r2 



/5a = 77^-^^ • (3-29) 



Notice that Pe = ^lA, as dictated by gauge invariance. 



4 Conclusions 

Constrained differential renormalization is a scheme of differential renormalization 
which does not introduce the ambiguities inherent to formal manipulations of singular 
expressions. At the moment it has only been developed at the one-loop level, although 
in principle it could be extended to higher orders. Such extension looks rather involved 
and is under study. In this paper we provide the techniques and explicit expressions that 
allow to perform any one-loop calculation in a four dimensional renormalizable theory 
(in the Feynman gauge) with this method. The whole procedure has been implemented 
in a computer package that performs all operations automatically in momentum space. 
This package is available and its description can be found in Ref. fiot] . Constrained dif- 
ferential renormalization respects the (one-loop) Ward identities of gauge invariance. 
Here, we have shown it for Scalar QED, which contains all the singular functions ap- 
pearing in renormalizable theories. In Ref. it was verified that the Ward identities in 
QED are satisfied and the correct triangular anomaly recovered (a more complete dis- 
cussion of anomalies will be presented elsewhere). That non-abelian gauge invariance 
is preserved is explicitly shown in Ref. |Q. The method also preserves supersymmetry 
in the calculation of the supergravity corrections to the 5^ — 2 of a charged lepton [0, |8| . 
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A Appendix 



In this appendix we derive the renormahzed expressions of all the singular basic func- 
tions appearing in renormalizable theories. Essentially, we use the techniques of point 
separation and point contraction, although in the simplest cases some intermediate 
steps can be avoided. For simplicity, we deal first with massless functions and then 
discuss the minor modifications of the process needed for the massive case. 



A.l Massless basic functions 

Let us start with the one-point basic functions. We saw in the text that rule p!a] implied 
A-^[l] = aOd{x), with a arbitrary in principle. Consider now the expression A^[l]5(y). 
Using rule 

A^[l](x)%) = [A{x)6{x)f5{y) = [A{x)S{x)5{x + y)f 

= A^[l]{x)5{x + y) . (A.l) 

Integration on x leads to the equation 

6iy) J d4xA^[l](x)=A^[l](y), (A.2) 

which, using rule ^, requires a = 0. Therefore, 

A^[1]=0. (A.3) 

Similarly, locality and power counting imply that 74^[5^] = [d^A{x)5{x)]^ is of the 
form 

A'^id^] = {a'a + ^,'^)d,,6{x) , (A.4) 
and the same argument used before implies a' = fi' = 0, and hence 

A^[a^]=0. (A.5) 

So, the massless one-point basic functions are zero in CDR. (The same results are 
obtained with the general form of point separation and point contraction.) 

Consider now the two-point basic functions. B^[l] was already given in the text. 
B^[9^] is easily obtained using the Leibnitz rule and the fact that the two propagators 
involved are identical: 

B^[a,J = [Aix)d^A{x)f = [A{x)A{x)f - [d^A{x)A{x)f (A.6) 

and then 

B'^[n] is directly obtained from A'^[l] using rule|^: 

B^p] = [A{x)aA{x)f = - [A{x)6{x)f = -A^[l] = . (A 



{A{x)ff = ld,B^[l] . (A.7) 
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This is the simplest example of point contraction. The renormalization of B[9^c?;^] is 
more involved. First, we notice that its most general renormalized expression is 

B^[5;.a.] = \{d^d, - \6^,u)B^[l] + -^[fd^^d, + + , (A.9) 

where the first term results from solving a differential equation at x 7^ and /, g and 
Ijl" are arbitrary coefficients with dimension 0, and 2, respectively, parametrizing the 
ambiguity in the local terms^ Second, applying rule ^ to the delta function in A [5^] 
we obtain 

= A«[a^] = -d^B^'p] + B^[9,,n] = B«[a^n] , (A.IO) 

and using the Leibnitz rule to move the derivatives to the first (identical) propagator 
and Eqs. (|A77|) and (0 



B«[5^n] = -ia^nB^[l] + lOpB^'id^d,] - B^[a^n] . (A.ll) 



Inserting Eq. ( A.IO ) into Eq. ( A.ll ), we find the consistency equation 

- ^n^^B^ll] + 25pB«[5^5p] = . (A. 12) 

If we substitute Eq. (|Al9| ) into Eq. (|Al^ ) we find /x" = and g = —f ■ For the moment, 
/ remains arbitrary. It will be fixed from the relation between B^[df^di,] and T^[djj_dvdp] 
later on. 

The renormalization of three-point basic functions containing one d'alambertian 
is directly obtained via point contraction, which relates them to two-point functions 
already renormalized: 

T«[n] = -B^[l]5(x + y), (A.13) 
T«[na^] = -dl{B^[l]{x)5{x + y)) - B«[a^](x)5(x + y) = 

-\{dl + dl){B^[mx + y)). (A.14) 

For T functions with more complex index structure we apply point separation to the 
corresponding B function with one less derivative. T^[dp,di,] can be decomposed into 
a traceless and a trace part, plus a possible ambiguous local term (symmetric under 
and X ^ y and with the correct dimension) 

T^[a^a,] = T[dpd, - + ^VT'^fn] + ■^b6p,6{x)5{y) . (A.15) 

The traceless part is finite and the trace has already been renormalized, so we just have 
to fix the dimensionless coefficient h. Separating points, 

B^[5^](x)<5(y) = -nVTld^] + 25^T^[9^ap] - T^[n5^] . (A.16) 



^ The dimensionful coefficient ji" is forbidden from ruie but we stiail instead use recurrence 
reiations and ttie initial condition A^[l] = to prove tfiat it vanisties, tfius sfiowing that rule ^ is 
consistent with the other ruies. 
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The simplest way to solve this equation for b is to integrate on x (see Ref. |^). Integra- 
tion on the "separated" point, y, leads to a tautology (this is a general fact). We obtain 
the value b = The same process can be applied to T^[d^d,ydp]. Its trace-traceless 
decomposition is 



X ( -B«[l](^(x + y) + ^c5{x)6iy) ) , 



(A.17) 



where we have used Eq. ( A.14|) . Point separation of B^[9^5i/], followed by point con- 
traction of the T functions with one d'alambertian, leads to 



(A.18) 



B^[a^5,](x)j(y) = -D^T^ia^a,] + 2a^T^[a^a,ap] 

+ {d;dy + dyd:){B^[l]{x)6{x + y)) 
+ B'^[dpd,]{x)6{x + y). 
Integrating on x after inserting Eqs. ( |A.15| ) and ( |A.17D , we get 

= -ny J d^xT[dpdu - + / d^xT[dpd,dp - ^(V^p + + Supdp)n] 

5iy), (A.19) 



+ 



(^ + ^-/)V°' + (^ + /W) 



which provides two equations (one for each tensor structure) that determine the coef- 
ficients / and c. Calculating the (finite) integrals we find f = ^ and c = — |. 

For four-point basic functions, which are at most logarithmically singular, the 
situation is similar. Q functions with at least one d'alambertian are directly contracted 
into T functions: 



-□^ (T[l]<5(x + y + z))- 2d; iT[dp]ix, y)5{x + y + z)) 

-T'^[n]{x,y)5{x + y + z), (A.20) 
-d^dt (T[l]5(x + y + z))-dl (T[a,](x, y)5{x + y + z)) 
-dl{T[dp]{x,y)5{x + y + z)) 

-T'^[d^d,]{x,y)5{x + y + z). (A.21) 



Finally, the renormalization of Q\dp^dudpd(j\ requires again the combined use of point 
separation and point contraction. The trace-traceless decomposition is not unique in 
this case. The simplest possibility is 



1 



+ -^{SpiySpa + SppSjya + ^paKp) 



167r2 



d6{x)6iy)6{z)) 



(A.22) 
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Point separation of T^[df^d,ydp] gives 

T^'id^d^dpjix, y)5{z) = -u^q[d^d,dp\ + 2diq''[dpd,dpd^] - q^'pdpd.dp] , (A.23) 

and point contraction of (^^[^dp^dydp] plus integration on x,y give the equation 

= -n^l d^xd^yq[dpd,dp\ + 2di j dWyQ^ [9^5,9^5,] 

+ j dST''[d^d,dp\{y,z). (A.24) 



Inserting Eq. ( |A.23 ) and the expression obtained for T^[dp,dydp], and calculating the 



integrals, we find d = |, what completes the renormalization of Q\dpdudpdcj\ 
A. 2 Massive basic functions 

The massive case is analogous to the massless one. The degree of the leading singular- 
ities is not changed by the inclusion of masses. Point separation and point contraction 
are used in the same manner, although now the resulting expressions contain additional 
terms proportional to the masses (see Eqs ( |2.15| ) and ( |2.16| )). In order to use (when 
possible) symmetrically the Leibnitz rule, we decompose the basic functions into two 
pieces, one symmetric in the masses — that is treated as in the massless case — and one 
antisymmetric — that is less singular. 

Let us illustrate the new features of the general massive case. The renormalization 
of Am[l] = ^m{x)5{x)^ is a bit more involved than that of A[l]. The renormalized 
expression must be of the form 

A^[l] = {an + f,^)5{x). (A.25) 

The same argument used for A^[l] forbids the term proportional to 06{x). On the 
other hand, using rule ^ and moving derivatives to the massive propagator we get 

A^^[l] = -[A^(x)aA(:r)]^ 

= -aB^^ll] + 2dpB^^[dp\ + A^[l] - m^Blll] . (A.26) 

Then, integrating this equality (with rule §), 



~2 2 
fi = —m 



j d'^B«,[l]. (A.27) 



So, we only have to integrate Eq. ( 2.13| ) and take the limit nii 0. (The integrals 



of expressions containing modified Bessel functions are easily performed with the tech- 
niques described in Appendix C of Ref. ^].) The result is given in Eq. pl^ ). Observe 
that fl is proportional to m, as it should for consistency with rule |l^ 

K^[dp\ is renormalized analogously to A^[9^]. B^^^^i^] is given in the text, 
whereas B^^^^l^^At] is obtained generalizing Eq. ( |A.7| ) with the decomposition men- 
tioned above: 

Bmim,2 [^m] 



2:i 



1 

+ 



647r4 



(A.28) 



On the other hand, Bmim2['-'] weh as T and Q functions are renormahzed following 
the massless case but including the terms proportional to the masses, e.g.: 

Krm2 P] = ^K^rm2 \^\ " ^t, [1] , (A.29) 



to be compared to Eq. ( [A.§| ). Finally, Bmim2 

\d^dy\ requires some work to solve the 
differential equation. Afterwards, the local terms are fixed analogously to the B[9^9^] 
case. The final result is given in Table 3. 

B Appendix 

In this appendix we give all the ingredients needed to obtain the momentum space 
expressions of Green functions renormahzed with CDR. Since the renormahzed func- 
tions are well-defined distributions, they admit a finite Fourier transform without any 
regulator. The Fourier transform of a distribution /(xi, . . . ^Xn)-, where Xi are four- 
dimensional variables, is 



f{pi,...,Pn) = J d^xi . . . dV e^'^^-P' . . . e'''"-P-f{xi,X2, . . . , x„) . (B.l) 
Due to rule ^, total derivatives directly yield 

d"^'^-ipj^. (B.2) 

Delta functions give rise to reduced Fourier transforms, with linear combinations of the 
original momenta as arguments. For instance, 

f{xi, Xn-l)S{Xn-l + Xn) ^ /(pi, • • • ,Pn-l " Pn) ■ (B.3) 

Therefore, one only needs the Fourier transforms of the renormahzed basic functions. 
We shall perform Fourier transforms in the shifted variables Zj. The Fourier transform 
of basic functions that are directly finite without renormalization is simply a convolution 
in momentum space ||T7|: 



Ftl.mJO]ipi,...,Pn-l) 

d^k d{k) 



(2^)4 [k^ + ml] [{k - pi)2 + mf] . . . [{k - Pn-iY + ml_^] 

^tlm,...m^.AO]{Pl,---.Vn-l). (B.4) 
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The integrals 1^""^ appear in usual one-loop calculations in momentum space and can 
be calculated with standard techniques. 

To Fourier transform the renormalized basic functions of Tables 2, 3 and 4 we only 
need (besides Eq. (|B.4| )) the following integrals: 

d^x e'"" = log M V , B.5 

X p + m. 

d^xe'^' PKoimx) = Svr^— ^— - , (B.7) 

J d'xe-r^K,(m,x)0'^^^^^^ = 2.2.'ltl^,^AO]{p,p) , (B.8) 
j d4^,...p ^i-^i("^i^) 0Ko(m2x) = 32^^l(^)^^„ja](0,p) . (B.9) 

The Fourier transforms involving modified Bessel functions are easily obtained using re- 
currence relations among them [|^, |l8[ . The recurrence relations include delta functions 
to make them valid at the origin. 

In Tables 5, 6 and 7 we collect the Fourier transforms of the basic functions in 
Tables 2, 3 and 4, respectively. These basic functions in (euclidean) momentum space 
can be directly used in momentum space calculations (see Ref. []lO|). 



/ 



A^[l] = 

A«[a^] = 

BR[i] = 1 log^ 



B^[5^] = -iiPA*B«[l] 
B^[a] = 



Table 5: Fourier transforms of massless one- and two-point renormalized basic func- 
tions. The small imaginary part in the logarithms here and in the next table, allows to 
analitically continue these formulae into the Minkowsky region (p^ < 0). 
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A^[l] = ^m2(l-logf^) 

B^.n.. [1] = 1^ { ^ log ^ + log + Cm,m, ip) } 

^tm, [d,] = -hp, [B^,m, [1] + 16^^ H - + log ^(mf + + (^^^i^) 

+ {ml - m'^)Crmm2(.P))] 
~ 3{P^^P^' — 4^tMuP )Bmim2[l] 

— |^?Ti| Irn2mimi [fc/ifc;/ — -gS^i/k'^]{p, p) + ?T1-2 Irnim2m2 [^/i^i/ ~ 'g^fj.uk'^jiPiP) j" 

where 

Cmimiip) = (log(p^ + ml + ml + A^/^ + k) - log(p2 + + - A^/^ _ ^g)^ ^ 

A = + (mi + m2)^) (p^ + (mi - m2)^) 



Table 6: Fourier transforms of massive one- and two-point renormalized basic functions. 
The integrals 1*^^^ are defined in the text. The momentum k is the integration variable. 
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TmimamsP] = "^3 ImUams [1] (Px, Py) - ^m^m^MPx - Py) 

'^'mim2»n3 [^M^f^p] ~ Imlm2m3 [^^i^V^'p ^ \{^fivkp + ^^ip^u + ^upkp,)k'^]{px,Py) 



[^'^/ii/T,^^^„^2,„j,^ [Odp\{px,Py) + ^fj,p'^mim2m3 [^'^u]{Pxj Py) 



+ 57li^ {^tiu{Pxp+Pyp) + Sf,p{Pxu+Pyu) + ^up{Px p + Py ii)) 

Qmim2m3m4 ['-"-'] ~ ~1T^\^^im2m3m4[k'^]{pxjPyjPz) + Imim2m3 [1] (Px — PziPy ~ Pz) 
+ '2'Pzp^mim2mz[kp]{jPx-Pz,Py-Pz) - T^imama [°] (Po; -Pz,Py-Pz) 

Qmim2m3m4 ['-'^At^!^] ~ ~'ITT'4^rn[m2rn3m4[k^k^][px, Py, Pz) -\- Pz pPzy'^raim2m-i[^\{Px ~ PziPy ~ Pz) 
Pz p,^mim2m3[ki/]{px ~ PziPy ~ Pz) ~\~ Pzv^mi_m2mi\k^{px ~ PziPy ~ Pz) 
- T^im2m3 [^/^^'^l (P^ -Pz,Py- Pz) 

Qmim2m3m4 [^M^i^^P^o-] ~ ^m\m2m:)nn[k^kykpkrj — ^{Sf^u^pa + ^pp^va + Spa^i/p)k^]{px,Py,Pz) 
+ + + V'^'^p) (Qmim2m3m4[°°](f'^'f'y'P2) + 9^) 



Table 7: Fourier transforms of three- and four-point renormalized basic functions. Px, 
Py and Pz are the conjugate momenta of the coordinate variables x, y and z. The 
momentum-space basic functions on the l.h.s. arc assumed to depend on Px,Py (for T 
functions) and Px,Py,Pz (for Q functions). The integrals I^"-* are defined in the text. 
The momentum k is the integration variable. 



27 



References 



[1] D.Z. Freedman, K. Johnson and J.I. Latorre, Nucl. Phys. B371 (1992) 353. 

[2] P.E. Haagensen, Mod. Phys. Lett. A7 (1992) 893; R. Munoz-Tapia, Phys. Lett. 
B295 (1992) 95; D.Z. Freedman, G. Grignani, K. Johnson and N. Rius, Ann. Phys. 
218 (1992) 75; RE. Haagensen and J.I. Latorre, Ann. Phys. 221 (1993) 77; C. 
Manuel, Int. J. Mod. Phys. A8 (1993) 3223; D.Z. Freedman, G. Lozano and N. 
Rius, Phys. Rev. D49 (1994) 1054; J. ComeUas, P.E. Haagensen and J.I. Latorre, 
Int. J. Mod. Phys. AlO (1995) 2819; M. Chaichian, W.F. Chen and H.C. Lee, 
Phys. Lett. B409 (1997) 325; V.A. Smirnov, Int. J. Mod. Phys. A12 (1997) 4241; 
D. Anselmi, D.Z. Freedman, M.T. Grisaru, A.A. Johansen, BRX-TH-420, CPTH- 



S.553.0897, HUTP-97/A037, MIT-CTP-2666, |hep-th/9708042 . 



[3] J.I. Latorre, C. Manuel and X. Vilasis-Cardona, Ann. Phys. 231 (1994) 149; G. 
Dunne and N. Rius, Phys. Lett. B293 (1992) 367; V.A. Smirnov, Theor. Math. 
Phys. 96 (1993) 974; V.A. Smirnov, Nucl. Phys. B427 (1994) 325; V.A. Smirnov, 
Z. Phys. C67 (1995) 531; V.A. Smirnov, Theor. Math. Phys. 108 (1997) 953; O. 



Schnetz, J. Math. Phys. 38 (1997) 738; D. Prange, DESY-97-211, |hep-th/9710225 



[4] D.Z. Freedman, K. Johnson, R. Muhoz-Tapia and X. Vilasis-Cardona, Nucl. Phys. 
B395 (1993) 454. 

[5] P.E. Haagensen and J.I. Latorre, Phys. Lett. B283 (1992) 293. 

[6] F. del Aguila, A. Culatti, R. Mufioz Tapia and M. Perez- Victoria, Phys. Lett. 
B419 (1998) 263. 

[7] F. del Aguila, A. Culatti, R. Muhoz Tapia and M. Perez- Victoria, Nucl. Phys. 
B504 (1997) 532. 

[8] F. del Aguila, A. Culatti, R. Muhoz Tapia and M. Perez- Victoria, International 
Workshop on Quantum Effects in MSSM, Universitat Autonoma de Barcelona, 
September 1997, |hep-ph/97H'47l . 



[9] M. Perez- Victoria, UG-FT- 

[10] T. Hahn and M. Perez- Victoria, UG-FT-87/98, KA-TP-7-1998. 

[11] J. Schwinger, Phys. Rev. 82 (1951) 918; C.S. Lam, Nuovo Cimento 38 (1965) 
1754. 

[12] J.H. Lowenstein, Phys. Rev. D4 (1971) 2281; J.H. Lowenstein, Commun. Math. 
Phys. 24 (1971) 1; Y.M.P. Lam, Phys. Rev. D6 (1972) 2145; Y.M.P. Lam, 
Phys. Rev. D7 (1973) 2943; T.E. Clark and J.H. Lowenstein, Nucl. Phys. B113 
(1976) 109. 

[13] O. Piguet, S.P. Sorella, Algebraic Renormalization: Perturbative 

Renormalization, Symmetries and Anomalies, Berlin (1995), Springer. 



28 



[14] G. 't Hooft and M. Veltman, Nucl. Phys. B44 (1972) 189; C.G. Bollini and J. 
Giambiagi, Nuovo Cim. 12 B (1972) 20; J.F. Ashmore, Nuovo Cim. Lett. 4 (1972) 
289; G.M. Cicuta and E. Montaldi, Nuovo Cim. Lett. 4 (1972) 329. 

[15] P. Breitenlohner and D. Maison, Commun. Math. Phys. 52 (1977) 11,39,55. 

[16] G. 't Hooft and M. Veltman, Diagrammar, CERN 73-9 (1973). 

[17] F. del Aguila and M. Perez- Victoria, Acta Phys. Polon. B28 (1997) 2279. 

[18] M. Abramowitz and LA. Stegun, Handbook of Mathematical Functions, New York 
(1972), Dover. 

[19] G.J. van Oldenborgh and J.A.M. Vermaseren, Z. Phys. C46 (1990) 425. 



29 



